A new interpretation of dielectric data in molecular glass formers 
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Literature dielectric data of glycerol, propylene carbonate and ortho-terphenyl (OTP) show that 
the measured dielectric relaxation is a decade faster than the Debye expectation, but still a decade 
slower than the breakdown of the shear modulus. From a comparison of time scales, the dielectric 
relaxation seems to be due to a process which relaxes not only the molecular orientation, but the 
entropy, the short-range order and the density as well. On the basis of this finding, we propose an 
alternative to the Gemant-DiMarzio-Bishop extension of the Debye picture. 

PACS numbers: 64.70.Pf, 77.22.Gm 
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Broadband dielectric spectroscopy has developed into 
the most important tool for the study of glass formers. 
It is able to cover the whole relevant frequency range, 
from /iHz to THz 1 . Therefore it would be very desir- 
able to understand the dielectric susceptibility in terms 
of physical processes. In particular, one would like to link 
the a-peak of the dielectric data to the disappearance of 
the shear modulus at long times, the essence of the flow 
process. 

Such a link is in principle provided by the Gemant- 
DiMarzio-Bishop^ extension (GDB extension) of De- 
bye's treatment. The extension considers the molecule as 
a small sphere with a hydrodynamic radius th immersed 
in the viscoelastic liquid. The medium is characterized 
by a frequency-dependent complex shear modulus 



G(lo) = G^g^Lo), 



(1) 



where Goo is the infinite-frequency shear modulus and 
g(uj) is a normalized complex function, increasing from 
zero to one as the frequency goes from zero to infinity. 

For molecules with a weak dipole moment like OTP, 
the GDB extension is 



1 



with 



Qow - n 2 1 + c r g(iv) 



^ttGqo r 3 H 
k B T ' 



(2) 



(3) 



Here e(w) is the complex frequency-dependent dielec- 
tric constant (with the conductivity contribution already 
subtracted), n is the refractive index, ei ow is the low- 
frequency limit of e and T is the temperature. One needs 
only the knowledge of the molecular radius ru . Then one 
can calculate e from measurable quantities. 

For strongly polar molecules like glycerol and propy- 
lene carbonate, one should take the difference between 
the external applied electric field and the internal field 
seen by the molecule into account^.. In Onsager's 
scheme^, extended to dynamics by Fatuzzo and Mason 6 
and reformulated by Niss and Jakobseni 
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This complex quadratic equation still allows to calculate 
e(ui) from G(u>), provided r# is known. 
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(b) propylene carbonate 
160 K 
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FIG. 1: Gemant-DiMarzio-Bishop (GDB) expectation (con- 
tinuous lines) compared to (a) dielectric datap- of glycerol at 
196 K (b) dielectric datap- of propylene carbonate at 160 K . 
The better fits (dashed lines) are explained in the text. 

The molecular radius can be determined from NMR 
field gradient diffusion data for glycerol^ (th — 0.16nm) 
and for propylene carbonate 9 (rjj = 0.26nm) via the 
Stokes-Einstein equation. For glycerol, there is a dy- 
namic shear modulus measuremenljiS in the temperature 
range close to the glass transition. Using the shift factors 
of this measurement, one can calculate G(u>) at the tem- 
perature 196 K of a dielectric measurement. At this 
temperature^, n 2 = 2.26. Fig. 1 (a) compares the GDB 
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expectation of eq. J3J with these values to measured 
data. The calculation underestimates the peak frequency 
in e" by an order of magnitude, a discrepancy which has 
been noted earlier—. 

The same holds for propylene carbonate. There is a 
shear modulus measurement at 159 K, close to the 
glass transition. Fig. 1 (b) compares the GDB prediction 
(with n 2 = 2.19 14 ) to dielectric datai 5 . at 160 K. Again, 
we find the peak in e" shifted by a decade. Obviously, 
the undercooled liquid finds a much faster way to relax 
the molecular orientation than the Dcbye mechanism of 
a molecular sphere rotating in a viscous liquid. 

As pointed out by Niss, Jakobsen and Olsen 16 , one 
does not even get a good fit if one adapts the molecular 
radius, because if one adapts the peak in e", the high- 
frequency e' gets much too high. 

In order to find out which mechanism might be re- 
sponsible for the decay of the molecular orientation, wc 
compare the dielectric relaxation times to the ones deter- 
mined by other techniques, in the spirit of earlier compar- 
isons by Ngai and RendelUi, Blochowicz et al^ and by 
Schroter and DonthA 3 «. We recalculate all data in terms 
of a Kohlrausch- Williams- Watts (KWW) decay in time 
exp(— {t/rKww) 13 ), either by using pragmatical recipes 
from the literatureA 9,20 or by refitting the data. To get 
rid of the strong temperature dependence of the viscos- 
ity, the resulting KWW relaxation time is divided by the 
Maxwell time T M axweii = v/Goo- 

The choice of a Kohlrausch or KWW-function is prac- 
tical for the following reasons: 

(i) it very often gives a good fit 

(ii) the inverse of the absorption peak frequency in 
lu is close to tkww^ so one compares peak frequencies, 
independent of the stretching parameter /3 

(iii) for a shear modulus following the Kohlrausch func- 
tion 



t~kww 

^Maxwell 



(5) 



Usually, (3 lies between 0.4 and 0.6, so the ratio should 
be between one third and two thirds; the Kohlrausch 
relaxation time should be a factor 1.5 to 3 shorter than 
the Maxwell time 

(iv) for a shear compliance measurement like the 
one2i in OTP, the steady-state compliance J° obeys the 
relation 2 ^ 
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so for a Kohlrausch functional! 

_ r(2//?)/3 

r((i//3)- 
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Since the compliance measurement supplies all three 
values J°, rj and Goo, one can determine (3 and 

TKWW /TMaxweii without Calculating TMaxweii- 



As usual, the measured viscosity r\ of our three sub- 
stances is fitted in terms of a combination of two Vogel- 
Fulcher-Tammann-Hesse laws 



log 77 = log 770i 



Bi 



T - T t 



(8) 



with i = 1 and i — 2, respectively. The first of these 
two is valid below a temperature T\ , the second above a 
temperature T 2 < T\. Between T 2 and T\, one takes a 
linear interpolation between the two to ensure continuity. 



substance 


glycerol 


propylene carbonate 


OTP 


log(j7oi /Pa s) 


-7.1 


-8.92 


-11.89 


Bi (K) 


1260 


667 


1461.2 


T01 (K) 


118 


122 


178.4 


7\ (K) 


283 


193 


310 


log(?7o2/Pa s) 


-5.45 


-3.91 


-4.24 


B 2 (K) 


780 


191 


245.9 


T 02 (K) 


153 


150 


241.72 


T 2 (K) 


283 


175 


275 


T 3 (K) 


187 


157 


243 


G x {T g ) (GPa) 


4.58 


2.5 


1.6 


a (K- 1 ) 


0.023 


0.007 


0.0057 


6 (K- 2 ) 


2. M0" 5 


2-10" 5 


1.4-10" 5 



TABLE I: Viscosity and shear modulus parameters. Refer- 
ences see text. 

The infinity frequency shear modulus Goo is 
parametrized in terms of a Taylor expansion around the 
glass temperature T g 

Goo = Goo(T fl )(l - a{T - T g ) + b(T - T g f). (9) 

The parameters of these two equations for the viscos- 
ity and the infinite frequency shear modulus are listed 
in Table I. For glycerol^ and propylene carbonate^, 
the Vogel-Fulcher parameters were taken from viscosity 
data fits in the literature. For OTP, we fitted our own 
parameters to the many viscosity measurements in the 
literaturo 21 i 24 i 25 i 26 i 27 . In glycerol, Goo was fitted to the 
Brillouin shear wave measurement of Scarponi et aU&. 
For propylene carbonate, there is no Brillouin shear wave 
measurement. Therefore, the infinite frequency shear 
modulus had to be taken from a longitudinal Brillouin 
measurement 2 ^, assuming Goo = cn/4. In OTP, there 
is a shear wave Brillouin scattering measurement^, but 
the Goo-values from this measurement extrapolate to zero 
already at 308 K. Therefore wc took Goo(T g ) from this 
measurement, but determined the parameters a and b 
of eq. © from the combined evaluation of light and x- 
ray Brillouin scattering of Monaco et a<2i, assuming the 
same temperature dependence for the infinite frequency 
longitudinal and shear moduli. 

We start the comparison for glycerol in Fig. 
2. The normalized dielectric- 1 ! 11 ? 32 KWW relax- 
ation times are compared to those from mechani- 
cal data 10 i 12 i 13 i 33 i 34 i 35 i 36 i 37 , from dynamic heat capacity 
measurements 3 ^ 3 ^, from NMR 4 ^, from PCS (photon cor- 
relation spectroscopy)^, from TG (transient grating) 37 
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and from neutron spin-echo measurements at the first 
sharp diffraction peak— Dynamic light scattering (DLS) 
data^ 3 - (not shown in Fig. 2) tend to lie between di- 
electrics and mechanics, but otherwise the figure corrob- 
orates the earlier conclusion of Schroter and Donthi 3 -, 
namely that there seems to be a grouping into the faster 
mechanical relaxation and a slower heat capacity, dielec- 
tric, NMR, PCS, TG and neutron spin-echo relaxation. 
The mechanical relaxation times follow the expectation 
of eq. JSJ) within reasonable error limits. The others 
tend to lie a factor of about ten higher. With chang- 
ing temperature, both time scales move together with a 
roughly constant separation. This shows that the misfit 
of the Debye result is temperature-independent, unlike 
the deviations from the Stokes-Einstein relation at lower 
temperatures^. We will come back to this point in the 
discussion of OTP. Here, let us first discuss what one sees 
in each technique. 
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FIG. 2: Kohlrausch- Williams- Watts relaxation times in glyc- 
erol, normalized to the Maxwell time as described in the text. 
Symbols: open squares sheas-i; open up triangles sheaj2~; 
open down triangles compression—; open circles longitudinal 
acoustic— ; open circle with plus shear and compression"; 
open circles with cross open diamonds longitudinal 

Brillouin— ; plusses dielectrio-; crosses dielectrio—; asterisks 
dielectrio—; full up triangles heat capacitji— ; full down tri- 
angles heat capacitji— ; full left triangles NMR—; full right 
triangles PCS—; full circles TG— ; full squares neutron spin- 
echoS. 



In the case of mechanical and dielectric data, there can 
be large differences in relaxation time between a modu- 
lus and the corresponding compliance. This difference is 
negligible if the relative change of the quantity in ques- 
tion is small, but here we deal with large relative changes. 
Therefore one has to check whether mechanical moduli 
and dielectric susceptibility are the correct choice. 

For the mechanical data, there is a good physical rea- 
son to choose the moduli rather than the compliances, 
because this choice leads to practically the same relax- 



ation time for the bulk and the shear modulus 33,44,45 . 
In fact, glycerol was one of the first cases in which this 
equality was demonstrated by the longitudinal and shear 
ultrasonic data of Piccirelli and Litovita^ 3 . (the circle with 
a plus in Fig. 2 at 255 K). In Fig. 2, it is again demon- 
strated at lower frequencies by the good agreement be- 
tween the shear measurementsi2*i&24 and the compres- 
sion measurement of Christensen and Olsen— If one goes 
over to compliances, this good agreement gets losli— 

Similarly, in the dielectric case one should take the 
dielectric constant rather than its inverse. Otherwise, the 
good agreement between NMR and dielectric constan^2ii£ 
(which is natural because both techniques sample the 
molecular orientation) would get completely lost. 

The transient grating (TG) experiment^ measures 
both the damping of longitudinal sound waves and a 
longer structural relaxation tim o 47 i 48 . One does not get 
the longitudinal sound wave relaxation time directly, but 
one can extract it from the temperature dependence of 
the damping. With a fitted (3 = 0.5, the damping of the 
sound waves translates into the two circles with crosses 
at 280 and 290 K in Fig. 2. In this case, one sees the 
splitting of time scales within experiments with a single 
sample and with the same temperature sensor. 

In the structural relaxation of the transient grating 
experiment, the heat of the phonon bath transforms into 
structural potential energy, thereby expanding the sam- 
ple. The relaxation time of this process is intimately 
related to the relaxation times of the heat capacity and 
of the density, which in turn are related to each other. 
The latter relation has been discussed earlier in sev- 
eral papersi2i 5 A 5 iii2i— Photon correlation spectroscopy 
(PCS) measures the refraction index fluctuations on the 
scale of the wavelength of the light, essentially density 
fluctuations. So it is not surprising to find the struc- 
tural TG relaxation times close to those of PCS and heat 
capacity. The neutron spin-echo measurements at the 
first sharp diffraction peak sample the decay of the short- 
range order of the molecular array. Again, it is not un- 
expected to find them close to those of the density and 
the entropy. What is surprising is to find the dielectric 
and the NMR times in the same group, because we are 
used to think of them as a single-molecule property^ and 
not as a collective process. 

The idea of two different time scales (or an initial and 
final stage of the same process) is further supported by 
the different shape of the relaxation functions for the two 
groups of Fig. 2. The mechanical data have a decidedly 
larger stretching (Pkww ~ 0.4. .0.5) than the heat ca- 
pacity, the dielectric constant and the neutron spin-echo 
signal (Pkww ~ 0.55.. 0.7). 

The same splitting of time scales, though not for so 
many different techniques as in the heavily studied case 
of glycerol, has been found for propylene carbonate^ and 
has been discussed in the framework of the mode coupling 
theory 54 . Note that this time scale splitting is not the 
two-stage scenario of the mode-coupling theory, because 
both time scales move together with the Maxwell time. 
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In fact, in ref™ the a-process of the theory was not at- 
tributed to the slower, but to the faster process. 



100 



10 



0.1 



NMR, 



propylene carbonate 



+ +++ + + + + + 



DLS 



□ mechanical O 



O Q> O O Q o < 



150 200 250 

temperature (K) 



300 



FIG. 3: Kohlrausch- Williams- Watts relaxation times in 
propylene carbonate, normalized to the Maxwell time as de- 
scribed in the text. Symbols: open square sheai— ; open di- 
amonds longitudinal Brillouin— ; plusses dielectria— ; full dia- 
mond NMR 9 ; open circles with cross DLS^. 

With the parameters in Table I, one can again re- 
late the measured Kohlrausch relaxation times to the 
Maxwell time. As in glycerol, mechanical sheai^ and 
Brillouini* data in Fig. 3 show a decade faster decay 
than NMR* and dielectric^ measurements, while DLS 
data 55 lie in between. 

Fig. 4 shows again a heavily studied case, OTP. 
Mechanical measurements include a shear compliance 
studySi, longitudinal ultrasonic data™, a transient grat- 
ing experiment™, a transverse Brillouin measurement™ 
and a thorough analysis of longitudinal Brillouin light 
and x-ray scattering™ Of these, the shear compliance, 
the ultrasonic and the longitudinal Brillouin data follow 
the Kohlrausch expectation of eq. JSJ, but the transverse 
Brillouin data and the longitudinal sound wave part of 
the transient grating results do not; they show a sud- 
den rise at about 290 K. The reason for this deviation 
is clearly revealed in the analysis of Monaco et ai™ At 
290 K, the Johari- Goldstein peak merges with the main 
a-relaxation. At such a point, our analysis in terms of a 
single Kohlrausch function is bound to fail. 

Otherwise, Fig. 4 corroborates the results in Fig. 
2. Again, the structural relaxation time from the TG 
experiment, heat capacity 5 ^, PCS™*2* and NMRffl lie 
close to the dielectric data™. The neutron data at the 
first sharp diffraction peak^ lie a bit lower, but do still 
clearly belong to the upper group. The dynamic light 



scattering point; 



,59.63 



do not lie between the two groups 



as in glycerol and propylene carbonate, but have higher 
relaxation times than all the other experiments. 

In OTP, there is a rather convincing explanation of 



the NMR data in terms of a single-molecule pictured, 
describing them in terms of rotational diffusion which 
follows the Debye-Stokes-Einstein equation 



D, 



k B T 
6irr]rH 



-rirD 



(10) 



where D tra ns is the translational diffusion constant of 
the molecule and D ro t is its rotational diffusion constant. 
For continuous rotational diffusion, the relaxation time 
for the Legendre polynomials is 



T~L,rot 



1 



L(L + l)D r 



(11) 



where L is the order of the Legendre polynomial. For 
the dielectric signal, L = 1, but the NMR measurements 
referenced so far are all two-pulse sequences for deuter- 
ated molecules, where L = 2. The OTP data™ are well 
described with a hydrodynamic radius Tn of 22 nm, close 
to the value rji = 23 nm found in NMR field gradient 
measurements^* at temperatures above 1.2 T g . The val- 
ues are smaller than the expected van-der-Waals radius 
of 37 nm, but still not too far away from it. 
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FIG. 4: Kohlrausch- Williams- Watts relaxation times in 
OTP, normalized to the Maxwell time as described in the 
text. Symbols: open squares sheas—; open circles longitu- 
dinal acoustic™; open up triangles TG—; open down tri- 
angles transverse Brillouin—; open diamonds longitudinal 
Brillouin--; plusses dielectric— ; full down triangles NMR™; 
full up triangles heat capacity—; full right triangles PCS—; 
full left triangles PCS—; full circles TG™; full squares neu- 
tron spin-echc™; crosses DL S 59 ' 63 . 

If one lowers the temperature, the rotational relaxation 
time follows the temperature dependence of the viscos- 
ity, while the translational diffusion begins to deviate to- 
wards higher values. The same decoupling between trans- 
lational and rotational motion has been found in photo- 
bleaching experiments** with guest molecules in OTP 
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and has been taken as evidence for dynamical hetero- 
geneity. In these experiments, one observes an increase 
of the relaxation times with increasing molecular diame- 
ter as expected, giving additional support to the single- 
molecule concept. 

But the GDB extension of this single-molecule picture 
to describe the relation between G(u>) and e(ui), eq. J2J, 
does not work. At the glass temperature of OTP with 
Gqo = 1-6 GPa, one calculates a c r of 80 from eq. (3). 
This implies that the peak in e"(oj) should be at a fac- 
tor of 80 lower frequency than the one in G"(u>), while 
the experiment shows only a factor of ten. Again, this 
discrepancy has been noted before^. In this case, one 
cannot blame the difference between external and inter- 
nal electric field, because the dipole moment of OTP is 
very small. 

In some cases, one even finds the peak in e"(uj) rather 
close to the one in G"(w). In decahydroisoquinolinc 
(DHIQ) at T g , they lie only a factor of 1.6 aparfei^, in- 
stead of the factor of about 100 that one expects. 

Also, the single-molecule picture fails to explain the 
striking coincidence between dielectric and NMR relax- 
ation times on the one hand and heat capacity, density 
and short range order relaxation times on the other. 

We will pursue an alternative explanation for the data 
in Figs. 2 to 4, namely that the flow or a-process be- 
gins at short times with a breakdown of the mechanical 
rigidity (the lower half of points in Figs. 2, 3 and 4). A 
decade in time later, there seems to be a final process 
which equilibrates everything, the density, the entropy 
and the short range order (the upper half of points in 
the three figures). This final process equilibrates also the 
molecular orientation, an order of magnitude earlier in 
time than expected on the basis of the Debye concept. In 
fact, a recent aging experiment^ shows that the dielec- 
tric relaxation time is indeed the final aging relaxation 
time also in a number of other molecular glass formers. 

A possible way to understand such a process is to 
postulate a configurational potential energy which has 
a small fraction of shear energy, able to decay within the 
mechanical relaxation time, while the large rest is merely 
feeding the shear energy. In a physical picture, one di- 
vides the potential energy of a given configuration into 
a long range shear part and everything else. This "ev- 
erything else" is supposed to be harmonic, decaying only 
via the shear energy channel. 

This is similar to the physical mechanism of the De- 
bye process, where the feeding energy is the energy of 
the electric dipole in the electric field. For the configura- 
tional energy, one replaces the electric dipole energy by 
the mechanical energy of an harmonic oscillator. As we 
will see, this change leads to a different equation; the two 
cases are similar, but not identical. 

To formulate the concept quantitatively, let us con- 
sider a mechanical model, a small spring r in series with 
a frequency-dependent spring g(tv) = G(w)/Goo. The 
compliance of the combination is the sum of the two com- 



pliances. Normalizing this sum, one gets 



$(w) = 



1 + r 



l + r/ S (w) 



(12) 



This is the Fourier transform of the decay function of the 
configurational energy according to the postulate above. 
$(w) is 1 in the high-frequency limit and zero in the low- 
frequency limit; it is a modulus function. 

We further postulate that the decay of the configura- 
tion involves a complete reorientation of the molecules, 
so it is mirrored in the dielectric signal. In the dielectric 
susceptibility, one expects to see 1 — $(o;). Then 



eM - n 2 _ 1 - g(u) 
Qow~n 2 l+g(w)/r' 



(13) 



This relation differs from the Gemant-DiMarzio-Bishop 
relation, eq. @, by the 1 — g(ui) in the numerator. 

We used eq. ()13ll to fit the dielectric data of glyceroli*ii 
and propylene carbonate^. g(u>) was obtained by first 
fitting the dynamical shear dataiSii^ at the glass tran- 
sition in terms of a KWW function and then shifting 
this function to the required temperature using the shift 
factors of the Maxwell time. In glycerol, we also took 
the slight change of the Kohlrausch-/? of the shear with 
temperature into account. 

To get a good fit, it turned out to be necessary to 
leave n as a free parameter and to allow for a slight dif- 
ference of the shift factor (remember that the shear and 
dielectric data stem from different laboratories). These 
deviations, however, remained small and of the order of 
the differences between the fit values for the two differ- 
ent dielectric glycerol experimentsiii. For the glycerol 
fit shown in Fig. 1 (a), the shift factor difference corre- 
sponded to a factor 0.7 in relaxation time. The n-value 
was 1.67 instead of 1.50. For the propylene carbonate fit 
in Fig. 1 (b), the shift factor difference corresponded to 
a factor of 0.6 in relaxation time, and n was 1.99 instead 
of 1.48. The shear energy fraction r was 0.19 in both 
fits. At higher temperatures, r tended to diminish, more 
strongly in propylene carbonate than in glycerol. 

In glycerol, the so-called "excess wing" at higher fre- 
quencies has been extensively discussed^SiSLS. From 
agin g 6 ^' 6 and pressure experiments^, one forms the im- 
pression that the excess wing is an unresolved secondary 
relaxation or Johari-Goldstein peak. Eq. 1|13|) does not 
add to this evidence, but demonstrates that one deals 
with both a dielectric and a mechanical excess wing. In 
fact, once the parameters are known, one can use the 
equation to calculate the expected shear response from 
dielectric data. It remains to be seen, however, how 
far one can trust an implicit assumption of eq. 113fl . 
namely that the elementary relaxation processes behind 
the breakdown of the shear modulus have the same ratio 
of the mechanical to the dielectric dipole moment over 
the whole frequency range. 

To conclude, we propose an alternative to the Debye 
concept, because the Debye scheme is unable to account 
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for dielectric data in glycerol, propylene carbonate and 
ortho-terphenyl. The alternative is based on the finding 
that the dielectric relaxation time in glycerol and OTP 
is close to the ones for the density, the entropy and the 
short range order. The scheme provides a good fit for 



broadband dielectric spectra of glycerol and propylene 
carbonate. 

We thank Peter Lunkenheimer and Ernst Rossler for 
supplying their dielectric data and for helpful discussions. 
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